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ABSTRACT 

An automatic elevator controller for the final phase of an 
Instrument Landing System approach is designed using optimization 
theory and the practicality of the controller investigated. The 
problem is discussed and the assumptions stated, Then a mathe- 
matical model for the aircraft and a desired flare-out approach 
path are derived. The aircraft and approach limitations are estab- 
lished and the model is — ne and the 
Parametric Expansion Method provide the optimal control from which 
the controller is designed. A computer program is developed to 
investigate the controller, The results are inconclusive and a 


recommendation for further study is made. 
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1. INTRODUCTION 

This study is an attempt to synthesize an automatic control 
mechanism for the terminal phase of an aircraft landing using 
optimization theory. The problem was originally presented by 
F. J. Ellert and C. W. Merrian [2]. Dynamic Programming [1] and 
the Parametric Expansion Method [2] are used to design a simple 
controller. 

Section two describes the problem and presents the assumptions 
made to simplify the initial study. Originally the investigation of 
Ellert and Merriam was to be continued but after several weeks of 
investigation, it was found that their equations for the aircraft 
represented only the changes about an equilibrium condition [3] 
rather than the true signals. Equations representing only the changes 
are derived in Appendix A and are used in section three to obtain 
different equations for the actual aircraft signals. 

In section four a more realistic desired approach is defined 
and equations for the desired values of the measured signals are 
derived. The performance limits and system constraints are dis- 
cussed in section five. Section six covers the testing of the 
aircraft model. 

Section seven develops the theory to obtain the optimal 
control in terms of the measured signals, aircraft constants, and funce 
tions of time: The necessary first-order, linear differential equation 
to obtain the time functions for the optimal control are derived 
in Appendix B. Section eight outlines the mechanization of a 
controller using the optimal control equation. Section nine 


discusses the computer programming and the investigation 
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2. PROBLEM STATEMENT 

Landing an aircraft in inclement weather when the pilot has 
no visual contact with the runway until the last few hundred feet 
is a function performed at many airports by an Instrument Landing 
System, commonly referred to as ILS, 

The aircraft follows an electronic beam which is set to a 
glide path of approximately three degrees and provides an azimuth 
reference to the runway. As the aircraft descends, the beam resol- 
ution becomes poor and sizable altitude errors can occur. The three 
degree glide path is normally followed to a point within half a mile 
from the approach end of the runway. The pilot then takes over 
visually to land the aircraft or is waved off to try again, 

At this time, a choice must be made by the pilot to follow 
ILS and hold the rate of descent he has or to make a flare-out 
approach to touch down on the runway with a lesser rate of descent, 
The flare-out approach puts less stress on the aircraft and is a 
more comfortable landing. 

A controller is needed to accomplish the flare-out approach so 
that if visibility is restricted, this kind of landing may still be 
made, The following assumptions are made to simplify these first 
studies: 

a. The aircraft is laterally aligned with the runway. Only 
errors in the vertical plane will be considered. 

b. At one-half mile from the end of the runway, the aircraft 
will be waved off if altitude, rate of descent, pitch angle, and 
pitch rate exceed specified limits. 


ec, Wind effects will not be considered, 
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3. SYSTEM MODEL DERIVATION 


The first step in the synthesis of an aircraft control system 
is the development of a mathematical model relating the control 
(in this case, elevator deflection) to some measurable response 
variables, It is essential that any model chosen be tested to 
insure that it follows all known responses of the real system. This 
will be discussed further in section six. 

Figure 1 defines the aircraft coordinates, angles, and positive 
direction of elevator deflection. Equations (3.1) and (3.2) have 
been derived to describe the motion of an aircraft [2] under the 
following restrictive conditions: 

a. The equations are linearized with the assumption that 
deviations from an equilibrium flight path are small. Only the 
changes about this equilibrium condition are represented by 
equations (3.1) and (3.2). 

b. The glide path angle, y, is small enough that the small 
angle approximations, siny=y and cos y= 1 can be used. This 
approximation is valid for the landing paths to be investigated. 

cc. It is assumed that the aircraft velocity, V, is held 
constant. 

Note that should a controller be found that is feasible from 
this basic synthesis, a more realistic model of the aircraft should 


be used for further studies. 
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Figure 1, Definition of Aircraft Coordinates, Angles, and Elevator 
Deflection. 
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(3.2) ah(s) 
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where AXs) = Laplace transform of the change in pitch rate, A6 (t). 
AS (8) = Laplace transform of the change in elevator 
deflection, Aé (t). 
Ah(s ) = Laplace transform of the change in vertical 


acceleration, A h(t). 
K.> T 5» We: — are constants with different values dependent on 
the type of aircraft. 
Equations (3.1) and (3.2) can be combined to provide a transfer 


function relating the change in vertical acceleration and elevator 


deflection, 
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K vu? 
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a7 4 2w 8 + we” 


(3.3) dh(s) = eS) 


which can be written as a fourth order linear differential equation 
1V eo@e 2 @o 
(3.4) Ah(t) + 2eW_ 4h(t) “ W, Ah(t) = K vw 89(t) 


relating the change in the control variable to changes in the 
response variables: change in altitude,4 h(t); change in rate of 
ascent, Ah(t) 3 change in vertical acceleration, Ah(t) ; and change 
in jerk, ah(t). 
All of these signals are not directly measurable. Signals 
that are readily measurable in the aircraft are desired. These are: 
a. Altitude, h(t), which can be measured from a barometric 
altimeter or radar altimeter. 
b. Rate of ascent, h(t), measured with a barometric rate meter. 
c. Pitch angle, 9(t). 
d. Pitch rate, a(t), both measurable from gyros. 
These are expressed in (3.5) as an equilibrium condition plus a 


change about that signal. 


(3.5) x, (t) = h(t) + Ah(t) 


x(t) = h(t) + Ah(t) 
xX, (t) = 9 ., (t) + Ad (t) 
x(t) = 6,(t) + ae (t) 


where the subscript e denotes the equilibrium condition. 

A set of first-order differential equations can be developed 
from (3.2) and (3.4) for the changes about equilibrium and result 
in equations of the form 


(3.6) 2(t) = Az(t) + Bu(t) 
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where z(t) is a 4xl vector and 
z(t) = [ah(t), ah(t), ao(t), ad(t)] 
A is a known 4x4 matrix of aircraft constants 
B is a known 4x1 distribution vector of the scalar control, 
u(t). 
The derivation of (3.6) is given in appendix A and represents 
the change of the measurable signals, x(t), about the equilibrium 


condition, Equations (3.5) can now be written 


(3.7) x, (t) = h(t) + 2, (t) 
x» (t) = h(t) + Z,(t) 
x4 (t) = a(t) + Z(t) 
x, (t) #8 Cae + z,(t) 
The equilibrium condition, x(t), is chosen to be a constant rate 
of descent, as variations about an ideal ILS approach are to be invest- 
igated. x, (t) becomes 


(3.8) h(t) h(O) + h(oyt 


h(t) = h(0) 


6 (te) 6(0) 


6 (t) 0 


Then (3.7) can be expressed as 
(3.9) x,(t) = h(O) + h(0)t + z(t) = x, (0) + x,(0)t + z,(t) 


y(t) se h(0) + ze) = x,(0) + 2z,(t) 


Hy 


x, (t) e(O) + Z(t) x, (0) + Z(t) 


x,(t) = 2,(t) 
where z, (0) = z, (0) = 2 (0) = 0 and z,,(0) = x, (0). 
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A set of linear differential equations describing x(t) is 

(3.10) x,(t) = h(0) + Z(t) 

k(t) = 2,(t) 

Ra(t) = 24(t) 

Gy 2” 22 0c) 
substituting (A.9) in (3.10) gives 
(3.11) x,(t) = h(0) + z(t) 

x, (t) = Ay 9% (t) + 5323 (t) 

Xg(t) = 2, (t) 

x(t) = a,o2,(t) + a,,z4(t) + a4,2,6t) + b,u(t) 
using (3.9), (3.11) ean be written. 
(3.12) x, (t) = X(t) 


Ko(t) = ayg[xy(t) - x,(0)] + ayy[x,(t) - x,(0)] 


q(t) x, (t) 


N 


x, (t) Ay o[X, (t) = x, (0) ] + a, 3[%3(t) - x4 (0) ] 

+ a, 1%, (0) + b,u(t) 

which are the desired linear equations as a function of the signals, 
x(t), and the control, u(t). The system model is also dependent on 


the equilibrium values of the descent rate and pitch angle because 


of the chosen equilibrium condition. 
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4. DESIRED FLARE-OUT APPROACH 

At time t = 0 with the aircraft some horizontal distance, D, 
from the end of the runway, it is desired to have the aircraft leave 
the three degree glideslope, increase its rate of descent, then flare- 
out to arrive over the end of the runway at time, t = th» at a lower 
altitude than if it had remained on the glideslope beam and with a 


lower rate of descent, From time, t = t,, until time of touchdown, 


1 


t = tes the rate of descent is held constant and all accelerations 
are zero. Figure 2 shows this desired path. An exponential altitude 


path over the time interval, t, < t< tes would be more realistic 


ih 


but the constant rate of descent was chosen for simplicity. 


hyo he —— — 







ALTITUDE 
($t) 


~ ww” GLIDESLOPE 


~ 
~ 
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~~ 


ais 


hy(t.) 


~ 


[+ D —- ++ 1000. —— 
END OF 


TOUCH DOWN 
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| 
t | 15 ~ te t s Ts 
TIME (S€C) 
Figure 2. Desired Flare-out Approach 


Note: All desired approach variables will be subscripted, d, to 
differentiate from the actual aircraft model signals. 
The following assumptions are made: 


a. Origin of the ILS glideslope beam and the desired touchdown 
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point coincide 1000 feet down the runway. 

b. Prior to t = 0, the aircraft has been following an ideal 
ILS approach. 

From the problem definition and the general desired approach 
specifications, the known conditions must be compiled to derive 
equations for the desired measurable signals. 

First the known conditions at times t = 0, t = ti» and t = t. 
will be stated. Using these boundary conditions, equations will 
then be derived for the time intervals 0<t < t, and t.st<t,. 

Prior to time t = 0, the aircraft is descending on an ideal ILS 
approach, Since all measurable signals are continuous, the following 
conditions exist at t = 0; 

h(0), h,(0), g,(0), o4(0), yg(0) are constants and 

h,(0), h,(0), 6,(0), a, (0); ¥g (9) are zero. 


The flight path specified over the interval t, <t< t. is 


1 


the same as before t = 0 (constant rate of descent). Continuity 


requires the same conditions at t, and t. with different constants. 


The desired flare-out approach path over 0 < t< t, in terms 
of altitude, h(t), and derivatives can be approximated by a 
seventh order polynomial since there are eight boundary conditions, 
These are 
a 5 6 7 


4° + at + act + alt 


four at t = O and four att=t 


(4.1) h(t) 


1 
h,(0) + h,(0)t +a 


ft 


° ° 3 4 5 6 
h(t) = h (0) + 4a,t + Sat + 6a-t + Jalt 
ee 2 3 4 5 
h(t) as 12a,t + 20a,t + 30a¢t + 42a.t 
ite) <=” 2ha,t + 6@a ecu ee t? + 210a.t" 
d 4 5 6 if 
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where a, = (20y(0) + 15y(t,)] _V_ - _35_[h,(0) - h,(t,)1 
(t,)”  (t,) 
1 1 
a, = [45y(0) + 39v(t))] -V_ + _84_ [h,(0) - hg(t,)] 
(t vi (t )? 
1 1 
ag'@ [36y(8) + S4v(ty)) tn = cTOe hae, hlea)) 
(t,) (t,) 
a, = [10y(0) + lOy(ty)] -V_ + _20_ [h,(0) - hy(t,)] 


6 7 
(t,) (t,) 


For the interval ty <te< tes the desired flight path is a 


constant rate of descent. The desired equations are: 


(4.2) ha(t) = hy(t,) + hy(t,)[t-ty] 


h(t) = hy(t,) 
h(t) = 0 
h(t) = 0 


The desired pitch angle, 8 Ct); can be expressed as its equili- 


brium value, 8 (0); plus changes about © ,(0): 
(4.3) 64(t) = 6, (0) + A@ qtt)- 
Taking the derivative gives 
(4.4) 6,(t) = A6,(t). 
From (3.2), 
(4.5) 46,(t) = Ts M(t) + 1 Ah A(t). 
d —— d —— d 
V V 
Integrating (4.5) 
(4.6) 40.(t) = Ts Ahi(t) + 14h,(t). 
d ——— ——~ od 
V V 
The desired values of (3.5) can now be written in terms of the 


derived equations 
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(4.7) 


x4 46) 


Xo qt) 


Xq y(t) 


xq) 


{ 


469 + h(0)t +a 


h : 6 
h qty) + h(t, (e-8y)5.t, < tac te 


469) + 4a 


O<'t <t 


O<t 


1 


4 


sty 


6 (0) oa 


a (4a,t” + 5a 


Vv 


3 


& 


4 + a t> + at + act 


t 


4 


t + 5a_t 


fee To t oe 


Ss. (i2a,t + 20a 


5 


5 


f 


2 


tt + Gat? + Jat’), O<t<t 


Ogft,) t,st<t, 


- 


a 


ls (12a, t + 20a 


(24a, t + 60a 


Zz 
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t2 


5 


3 


6 


+ 6a t? a /aat 


3 


be) 
5 


+ 120a, 


t? + 30a 


? 


6 


6 7 


6 
7 3 


+ 210a, t ws 


_* + 42a, t yy 


2 


t + 30a," + 42art = 
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5. PERFORMANCE LIMITS AND SYSTEM CONSTRAINTS 

Aircraft constants, numerical values of the desired flight 
path, and variable limits must be specified before testing the 
system model, The T-28 aircraft was chosen for this investigation. 
Values of the aircraft constants in (3.1) and (3.2) for the T-28 


in a landing configuration are; 


(5.1) K, -1.350 enim 


T = 2.205 sec 

W. = 0.6161 radians/sec 

—E = 0.4230 

V = 157.8 ft./sec 

The constants for equations (4.7) are calculated from (5.1) 

with reference to figure 2. ILS weather minimums of visibility 
vary from a half to a quarter mile. Using this as a basis, a value 
of one-half mile is chosen for D. In consistent units, 
(5.2) D = 2640 ft. 
Then from figure 2, 


(5.3) ha(0) = (D+ 1000)tan 3° = 190.8 ft, 


Using the small angle assumption, 
(5.4) h,(0) = y,(0)V 


-0.05236 radians 


where ¥ 469) 


and 


(5.5) ty BD 
V 


t = (D f 1000) sec 
£ V 
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A value of twenty feet is chosen for the desired altitude 
over the end of the runway, then 


(5.6) h(t,) = 20 ft. 


Yqg6ty) = «0,02 radians 


hy (t,) = ght pv = «3,156 ft./sec 


The physical limitations of the aircraft must be considered; 
important details are discussed below. 

1) Angle of Attack, The T-28 aircraft stalls at an angle of 
attack, a(t), of approximately 21° and a speed of 7/2 knots. As 
the angle of attack approaches this limit, the speed cannot be main- 
tained constant and the linearized system equations no longer repre- 
sent the aircraft motion. It is assumed that realistic representation 
is lost for a(t) > 138°, This value will be called the stall angle of 
attack, ae For the T-28 in landing configuration descending on an 
ILS glideslope at 93.5 knots, the angle of attack, a(0), is 11° or 
0.1920 radians. 

2) Pitch Angle. The pitch angle at touchdown must be between 
the limits 0° < 6(t .) < 14°, The lower limit prevents the nose 
wheel from touching first and the upper limit prevents the tail skag 
from dragging on the runway. 

3) Elevator Control. The elevator motion is limited by the 
mechanical stops of the actuator arm. Assumed values are -35° and 
+15°, The direction is defined in figure 1. Aircraft motion with 
the elevator against these stops for non-zero time is not permitted 
for the linearized model used, 

From figure 1, 


(5.7) O(t) = a(t) + y(t) 
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and the desired pitch angle at t = 0 is 
(5.8) 8, (0) - ag (0) of yg 69) = 0,1396 radians. 

Since the angle of attack must remain below the stall value, 
a, an equation for a(t) is desired, From (3.2), the change in 
pitch rate is 


(5.9)  ao(t) = Ts ah(t) + _1. sm(t). 
V V 


(5.7) may be written in the form 


(5.10) 6 (t) + ao(t) = a(t) + dat) + y(t) + ay(t). 


Taking the time derivative with 0, (t) = a(t) = ett) = QO, 
(5.11) ae(t) = A&(t) + ay(t) 
and equating (5.9) and (5.11) gives 


(5.12) Aa(t) + Ay(t) = _Ts ah(t) + _1 ah(t). 
V V 


It can be shown [3] that Aa(t) = TAY so the T, term of (5.12) is 
associated with Aa(t) and 


(5.13) A(t) = Ts Ah(t) 
V 


(5.14) Ay(t) = _1 Ah(t). 
v 


Integrating equations (5.13) and (5.14) yields 
(15) a(t) c= o Gt) + Ts Ah(t) 
(5.16) y(t) = y(t) + _L h(t) 
V 
where a(t) = a(0) and y ,(t) = y;(0). 


From (5.4), y(0) # h(0) and (5.16) becomes 
V 
(5.17) y(t) = _1 x)(t) 
V 


Using (5.7), %(t) can now be determined from the system model 
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by 


(6:18) a(t) = x, (t) - a x, (t) 


In section 2, it was stated that if the measurable signals, 
x(o), were outside specified limits, the aircraft would be waved 
off for another approach. These limits are somewhat arbitrary and 
should be chosen only after a close study of existing conditions of 
actual aircraft making ILS approaches. 


For this study the following limits are assumed: 


(5.19) x, (0) = x1 60) + 20, ft. 
x4 (0) = Xo q (9) + 2 ft./sec 
x, (0) = *, 00) + 0.0350 radians 
x, (0) = x, 4 (0) + 0.0043 radians 
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6. MODEL TESTING 

Every mathematical model derived for a physical system should 
be tested by applying known controls to insure that the model 
accurately describes the system, Originally the z(t) equations were 
used as the aircraft model assuming that they represented the 
signals, x(t). In this problem there are two controls that can be 
used for testing. The control, u(t) = 0, applied to the model allows 
the aircraft to follow its equilibrium flight path which is a constant 
rate of descent to the runway. From (3.4) a control can be derived 
using the equations, x,(t), since they are continuous and known 


functions of time. 


VV ee0o 
(el) a(t) me kit) + _2e Ah(t) + a he 
K_VW ‘ AL al 
s§ 8s 
Since h(t) and h(t) equal zero, 
(6.2) m(t) = y(t) 
Ah,(t) = h(t) 
Using (4.1) and (4.2), 
(6.3) h(t) =" I2e,t° + 208.0°> + Jomt* + wear” 
i d 4 5 6 7 
Ah,(t) = 24a,t + 60a ce + 120a te + 210a e 
d 4 5 6 7 
IV 2 3 
A hy (t) = 24a), le 120a,t + 360a,t + 840a,t 


forO < t < ti. 


as ooo IV 
(6.4) Ah y(t) = dh ,(t) = dh ,(t) = 0 


fort. < t << t 


1 Ee 


With (6.3) and (6.4) used in (6.1) over their respective time intervals 
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and this u(t) applied to the model, it should follow the desired 
flight path shown in figure 2. 

When these controls were applied to the first model, it did 
not follow either flight path (see figure 3). After deriving the 
model in section 3, the same two controls were applied and the 


proper flight paths were followed (see figure 4). 
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7. APPLICATION OF OPTIMIZATION THEORY 

The first step in applying Dynamic Programming is to obtain 
an error index, This index determines the form of the resulting 
optimal control. The errors, x(a) - x,(4), for Ounm @ < te; 
the errors at the final time, x(t.) ~ x, (t_); and the control, 
u(t), are considered to be significant for this problem. The Bolza 


form of error index is considered adequate to insure satisfactory 


System performance. The error index, J, is defined by 


GD J = 1 fixltp - alto)" Fixed - x4(t91} 
t 
+ if * Ftc) = xy(ed)7Q¢e) [xed ~ x4(t)] + u?(e) } at 
Oo 


where F is a 4x4, positive, semi-definite, constant matrix. 
Q(t) is a 4x4, positive, semi-definite, time variant matrix. 
Both F and Q(t) are arbitrary weighting functions of the 
System errors and are chosen to be symmetric. 


The condition by which (7.1) is a minimum is given by [1] 


(7.2) min F { [x(a) = x,(a)] Qa) (x(a) - x3(a)] + "(ap ae| = 0 
u(a){2 da 


where a is an arbitrary instant of time in the interval,O < a < te. 


ty is the fixed final time. E is defined as the minimum error 


function, E{f(a), x(a)] = min J, which is assumed to be only a 
u(a) 


function of a and the signals measured at a, 


As a=~t it can be shown that 


f? 
(7.3) E[E(t gd» HCE] = Lfix(tp) - xg(t_)] Fix(t,) - x4(t.)]} 


Define 


(7.4) tx(a), w(ad] = LffeeCa) ~ xy(2)]7QC@)[x(@) - xy(@)] + uv? ¢ay} 
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then (7.2) can be written as 


(7.5) min [HIx(@), aL + au = 0 
u(a) da 


The derivative of E with respect to a may be written as 


(7.6) dE = d Ef[f(a), x(a 
da da 
4 
= DE df(a)_ + >, 9Xn QE[f(a), x(a)] 
@f(a) da n=1 da dx, (a) 


The minimal error condition (7.2) can then be written 


(7.7) min [Ftx(a), u(a)] + YE ___d£(a) + (ORE 6 
u(a) “Sf(a)_ da ax, 


x(a) XE. + X,(a) DE + x,(a) XE | = 0 
, aX, . 355 ° | 


For (7.7) to be a minimum for all a, 0 < a < t,, the first 


f? 
derivative of the term in brackets with respect to u(a) must be zero 
and the second derivative must be positive. E was defined as a 
function of a and x(a) and is not a function of u(a). From the 
system equations (3.12) it is seen that x x, (a) is the only time 


derivative that is a function of u(a). The derivatives of (7.7) 


with respect to u(a) can then be written 


(7.8) — » *4(a) 
ah ~ yuca) 


NH +p JE X46 (a) > O 
aee)- 9) du 
From (3.12) and (7.4) 


(7.9) o%4(a) = by» 4°%4 (a) =, © 
‘bufay Mula) 
)_H = u(a), »H = 1 
Yu(a) \ ula)? 


Substitution of (7.9) in (7.8) yields 


(7.10) u(a) = ~b, bE 
d 


1 > 0 
(7.10) defines the optimal control for this system over the interval, 


O<ae< tee Substituting (7.4), (3.12), and (7.10) in (7.7) gives 


2 
(7.11) 1 ftx(@) - x,(2)]" Q(a)[x(a) - x,(a)]}+ ae + 
2 


9 


df(a) + x,(a) E + (a) - x,(0)] + 
3) da 3x, [222l%9(8 age 


45313 (4) “ 3(0)]| QE + x, (a) QE + feuatxa6e) ~ x» (0)] 
dx, dX, 


2 
+ a,3[%3 (a) © x3 (0)] + ayiy(2)| QE ~ aed = Q 
4 4 


which is a partial differential equation defining E. If (7.11) can 
be solved for E, then the optimal control of (7.10) can be found. 

The Parametric Expansion Method [2] is now applied to solve 
(7.11) and provide an optimal solution. The method is outlined as 
follows: 


Assume a form for E consistent with the functional description. 


4 4 4 
(7.12) E = k(a) - 2 = k (a) x(a) +2 2 Kap (% (2) p (2) 


where k_ (a) = k_ (a). 

mp?) = kn(@) 
Taking the derivative with respect to x, (a), substitution in (7.10) 
yields the optimal control as a function of time, the aircraft para- 
meters, and the measurable signals, x(a). 


(7.13) u(t) = 2b, [k, (t) - k(t), (2) - ko, (t) x, (t) 


- ky, (t) x, (t) - k,(t)x, (t)] 
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since a is an arbitrary instant, it can be real time, t. Obtain 
the derivatives required by (7.11) from (7.12). Substitute these 
in (7.11) and collect terms so that the expression is in the form: 
(7.14) PB, + Pyx,(a) + Pox, (a) + Pyx4(a) + P,x,(a) + Py x, (a) 

ir Pox 9 gk + P33%3 te + PX, (a) 4 P yo (4) Xo (a) + 

P 3% 1 (4)%3 (a) + Pi 4%) (2)*, (4) + Po 4Xo(a)%4(a) “} 

Pou x, (a) x, (a) + P34% 3 (ax, (a) = O 
The P coefficients contain terms of the k-parameters, time 

derivatives of the k(a)'s, the desired signals, the weighting factors, 
the aircraft constants and the chosen equilibrium conditions. For 
(7.14) to be valid for all values of the measured signals, each P 
coefficient must independently equal zero. 


(7215) Po os Ge = 0 
fe) m mp 


where m = 1, 2, 3, 4; p = 1, 2, 3, 4. 
Equation (7.15) gives a set of 15 independent, ordinary, first- 
order differential equations which define the k(a)'s. These equations 


with the boundary conditions are derived in Appendix B. 
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8. IMPLEMENTATION OF A CONTROLLER 

Several types of controllers could be designed depending on 
the assumed form of E and the methods available for measuring time 
and velocity. The block diagram shown in figure 5 is obtained from 
(7.13) and is a simple analog control system, A constant speed motor 
rotates the shaft at a constant velocity. The k(t)'s of (7.13) 
represent the settings of wire-wound potentiometers where the wiper 
arms rotate with the shaft. The measured signals are multiplied 

by the appropriate k(t), summed and multiplied to provide the ele- 


‘vator deflection signal. 


e 


lov h (t) h (t) 6 tt) 6 ¢t) 


CONSTANT | © 











EACH POT IS WIRE-WOUND 2 
AS A FUNCTION OF THE S(t) TK,W, 
RESPECTIVE k(t). 


Figure 5. Elevator Controller 
In the derivation of the k(t) equations, the equilibrium values 
of the descent rate, Xo (t) » and pitch angle, Xq(t), must be known. 


This means that the optimal k(t) equations are dependent on the 
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equilibrium conditions, This is a function of the aircraft model 
and would be altered with a different aircraft model. Little error 
should be introduced if these values are assumed to be the desired 


descent rate and pitch angle at problem start. 


(8.1) X(t) Xo 4 (0) 


x, (t) x, 4 (0) 


By assuming these values the system is no longer optimal, but it is 
desired to sacrifice optimality in the interest of simplicity and to 
investigate the feasibility of this simple controller. The k(t) 
equations are also dependent on the aircraft constants so a diff- 
erent set must be obtained for each type of aircraft. 

At this point two big questions remain. 1) Can k(t), k,,(%)> 
kay ft) k,,(t), and kyy ft) be found such that, for all aircraft 
with x(0) within the signal limits given by (5.19), the aircraft 
will make a satisfactory landing? 2) Are the k(t)'s found above 
linear enough to be wound on potentiometers? A method to investi- 


gate the answers to these questions is outlined in section 9. 
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9. INVESTIGATION PROCEDURE 

Since all of the constants and variables of the k(t) equations, 
(B.5), have been calculated or are known except the weighting factors 
F and Q(t), the answers to the questions in section 8 can be investi- 
gated by the following procedure: 

1) Arbitrarily pick the weighting factors of the error index, 
F and Q(t). 

2) Integrate equations (B.5) from t = te to t = 0 to obtain 
the k(t)'s of (7.13). 

3) Integrate the aircraft model equations (3.12) from 
t= QO tots= t. using u(t) of (7.13) with the k(t)'s obtained in 2. 

4) Observe the final signal errors, x(t,) - x,(t,). If these 
are not within prescribed limits, change F and/or Q(t) and repeat 
steps 2 through 4. 

A program was written to accomplish the above procedure using the 
CDC 1604 computer and an existing modified Runge-Kutta integration 
subroutine, There were two main sections in the computer program. 
1) Equations (B.5) were integrated from t = t, tot= 0 with the 
values of k(t), k,(%)> k(t) » k(t)» and ki Ct) stored in a 
buffer. 2) Three sets of equations (3.12) were integrated from 


t = 0 to t = t. using u(t) of (7.13) and the k(t)'s that were 


f 
stored in the buffer. 

Three aircraft models were chosen for each program so that one 
could be started with positive errors, one started on the desired 


trajectory, and the third started with negative errors. This allowed 


all of the system errors to be observed in one program. 
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The final error limits were arbitrarily chosen as: 

a) An altitude, descent rate combination such that all three 
aircraft land within +500 feet of the desired touchdown point with 
a descent rate of + 100 ft./min or + 1.66 ft./sec from the desired 
rate of descent, Xo q(t_)- 

b) Pitch angle as stated in section 5, 0° < Q(t -) < 14°, 

c) No limit was set on the pitch rate but it must be small 
for the approach to closely follow the desired flight path. 

For this investigation, the aircraft was considered to have 
only an altitude error and Q(t) = Q, a constant matrix consisting 
of the diagonal elements only. The study was conducted in the 
following manner: 

1) Determine if successful approaches can be made with con- 
Stant weighting of any one signal; or any two signals; or combin- 
ations of any three signals; and finally, combinations of all 
four signals. 

2) Determine if successful approaches can be made with constant 
weighting of all signals, progressively adding diagonal elements of 
the F matrix. 

3) When successful approaches are found, determine the nature 


of the k(t)'s required for the controller. 
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10. RESULTS 

Digital simulation was used to study the areas outlined in 
section 9. The results were: 

1) Constant weighting Of all four signals was required to 
obtain a satisfactory trajectory. All final error limits were 
met except that of the minimum descent rate. Figure 6 shows the 


altitude vs. time/distance for three aircraft started with an 


initial altitude error of + 10 feet, where 17 oc at Mowe 


0 4 


“= 8Jsiiol_-, dg, = 6.0x10"°, and a, “ 88x10". All 
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other q's and F were zero. The minimum descent rate at touchdown 
was -5.36 ft./sec compared with the limit of -4.91 ft./sec. 

2) This area was not fully investigated but successful 
approaches were made with constant weighting of the four signals 
(Q equal to zero except 4112 400° 933: qdy4 and weighting all 


four signals at the final time (F equal to zero except f ft 


bl! 22> 


£ Figure 7 shows the altitude vs. time/distance for three 


33° f4y)° 
aircraft started with an initial altitude error of + 20 feet where 


-10 =5 
P | 9992 9 933 Se 7.00x10 9 
7 


aa a -9 Ss - 
8.00x10 " and £,, = 1.80x10", f,, 2.00x10 ‘, 
2,.90x107>, £) = 3,00x10°>, All final errors were within 


7.10x10 = 2ee5enD’ 


rh 
i 


specified limits. 

3) For the approach of part 1, where F is zero, k, (t) is shown 
in figure 8, k(t) in figure 9, and Kk, 1,6) in figure 10. kt) and 
k,,,(t) were insignificant.(i.e. nearly zero). 

For the approach of part 2, k(t) is shown in figure 11, k(t) 
in figure 12, and k, 1,6) in figure 13. Again k 46%) and ko, Ct) were 


insignificant. 


39 


(208) owpZ 





si ie SA 2b 


°9 qaHNoOld 


“HONI7SLINM G@+3OOT = AIWIE-A 
"HONIG LINM BB+IVO'C = JAIWIE-% 


T60- 


(az) epnataty 


40 


a TQ 
bie ice 


i Se 


‘£ WANSIA 


"HONIZELINA CO+3VOT = JWIS-1 
ONT /E LIN BBFAVO"C = IP WIS-y 


(998s) out] 





(33) @pnazIaty 


41 


(99s) ewtL 


amo dua INTL “SA Ay 


°8 aanola 


"HONI/“S LIN TO-FBO'C = AWIC-4 
"HONI“S LINN BB+AOO"C = JWIC-»% 





KINO d 4d JIT | “SCY ACY 


°6 Fanos 


"“HONIYGLINT FQ-380'C = IWIC-h 
HONIYG LINN @B+300'C = IWIC-x 


(99S) out] 





(098) aut] 





Se ae et) n> 


"OT FANOLA 


“HONIZGLINA TQ-FeB'C = IWIC-h 
"HONI“G LIN BO+IOO"C = IWIC-X 


weeny AT) SO). Te 


‘Tr manora 


HONT/SLINT 1Q-308"C = JWOC-4 
HONI “SLING O8+380"°C = FWOC-%X 


(9968) sutl 





(99s) sutL 





Paw emt TL “SA ACD 


“Cl waNnsid 


"HONT/SLINA TQ-RWO'C = IWIC-4 
"HONT/SLING BB+300'C = JWOS-X 


v 


46 


(das) outL 





vow pe ATT SA nny 


"€T adnold 


“HONT “SLING TQ-308'C = JWIC-+ 
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11. CONCLUSIONS 

For the results specified in section 10, the following conclu- 
sions are stated: 

1) Although the minimum descent rate error was not within 
limits using the Q matrix only, it is felt that with several small 
adjustments of the weighting factors the descent rate could be brought 
within limits. The k(t)'s are of the same sign with maximum values 
less than one and could be wound on the potentiometers. 

2) While all of the error limits were met for the approach 
using both F and Q, the k(t)'s changed sign and are unsatisfactory 
for this controller. As adjustments were made to F and Q, bringing 
the flight paths closer to the desired trajectory, the k(t)'s 
approached the same sign over the interval of flight. It is believed 
that if F and Q are adjusted for flight paths with smaller errors, 
this problem would be eliminated, 

Two recommendations are made: 

1) That investigation be continued until definite conclusions 
can be obtained about the practicality of this controller. 

2) That a method be investigated to systemmatically adjust 


the weighting factors, the elements of the Q and F matrices, 
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APPENDIX A. DERIVATION OF z(t) EQUATIONS 
Define 


(A.1) z(t) = Ah(t) 


ii 


z(t) = Ah(t) 
z,(t) = AQ(t) 
z(t) = AO(t) 
u(t) = AS(t) 


The time derivatives of the z(t)'s are: 


(A.2) 2, (t) = MOh(t) Zo (t) 


z(t) = Ah(t) 
Z,(t) = Aect) = 2, (t) 
z(t) = A@(t) 
Equations for Ah(t) and A8(t) in terms of the states, z(t), are 


needed. (3.2) can be written as 


(A.3) T Oh(s) + Ah{s) = VAQ(s) 


Taking the inverse Laplace transform and solving for Ah(t) yields 


(A.4) Ah(t) = -1,ah(t) + _vAQ(t) 
T 


s s 
Differentiating (A.4) gives 


(a.5) Ant) = - 1Ah(t) + _VA0(t) 


wy, 
T 
8 s 
Integrating (A.4) and assuming Ah(0) = 0 gives 
(A.6) Ah(t) = -1Ah(t) + _VAQ(t) 
8 Ms 
Substituting equations (A.4), (A.5) and (A.6) in (3.2) and solving 


for AQ(t) gives 


30 


(A.7) AQ(t) = 1 
Vee: 
T Ss 


-fa - 28Ws + w.*] aocey + [_1 - 28_] adce) 
ae a ie ‘ 


2 
kel AS(t) 


le 2 “ An(t) 
Ss 


Inserting (A.6) and (A.7), written in terms of (A.1), into (A.2) 


yields 
(A.8) z(t) = 2,(t) 
z(t) = - 1 2,(t) + _V z(t) 
2 - 2 7 3 
Ss Ss 
2,(t) = z,(t) | 
z(t) = ee Siig + W. | zo(t) -|1 -_2BWe + WB (t) 
T Ts T 5 T 
$s & & 
+ a - z(t) + K.T.W, u(t) 
T 
Ss 


which is in the form 


(A.9)  2(t) = Az(t) +Bu(t) 





where 
0 1 0 =O 0 
Az |°% 402 493 9 v= 0 
0 0 oO 1 0 
O ano 443 844 oa 
#2 
and a). = -_1, 4, = Wy a. * 1] 1 - 28Ws + ws 
T T V 2 
s s a ly T 
s S 
a = =| 1 = 2§W stile a a 1 - 2$w bak TW 
43 a — s} ’ 4 ‘s's's 
T 
tT. T; S 
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APPENDIX B, DERIVATION OF 


ae EQUATIONS 


The partial derivatives of E defined in (7.12) are needed for 


substitution in (7.11). § 
the interval 0 < a <¢ t 
functional notation, (t), 


simplicity. 


The required partial 





ince a is an arbitrary instant of time in 


it can be chosen as real time, t. The 


£? 
on all the variables will be dropped for 


derivatives of E are: 


(B.l) dE df(t) = YE Ek) + JE k, 
F(t) dt k ak, dk, 
@ @ td @ e 
+h Ks + oe 47 ae Ku + 2 “22 ie 33 
3 4 ial 22 33 
) 6 & & é 
+ LE ki, + ut ko + tf ki, + SE ki, + bE kog 
44 12 13 14 23 
e 
+ ae ko, = Ka, 
24 34 
+ = 1, VE = 2X1» de = “2X05 Ae = -2%45 
My 2 3 
SE = 86 AR am ys a8 aE a me, 
Vk + Ye 1% 2. ak : 
4 ll 22 33 
2 
ue. = 2 . S = 2x), YE = 2x,%4; _ 2x,%) > 
44 12 13 14 
te ss 2x.%3+ dE = 2942 VE = 2x,X,; 
23 24 34 
XN OF 2(-k, + ky y%) + kyoX, + ky 4X3 + kyA%,)> 
Ue, 
= 2(-k, + k,4x, + k,ox, + ky 4X4 + k,,%,)> 
2 
2. = 2(-k, + ky 4x, + ky4x, + kggx, + ky,%),)» 
3 


ys 


sE = 2(-k, + ky 4%, + ky), + k3,%4 + k,,%,); 


2 
2 
fA 4(k, 0 = 2k, ky), 7 2k, ko, X, -2K, ky, x, 


Zz 2 Zn 2 2 2 2 
- 2k ky ,%, + Kia xy + kon Xo + ky, X4 + Kyy * 


F 2k ky % 1% + 2ky Akg, %1%3 + 2ky kK, x 1%, 


+ 2ko kg,¥o%y + 2ko,k,%oX, + 2ka,k, 4%, ) 


The derivatives of (B.1) are then substituted in (7.11) to obtain 
T 25 mee 
(B.2) Lye - X,) Q - z.)} “2b, (k,” = 2k,k, x, ~2k, ky, X5 


2 2 Pe? Zae 
2k 1k ,%3 2k kak, -} Ky xy + Koy Xo + kK, x4 


2 2 
Kay %y + 2KyyKog¥y% + 2K yk o,%1%3 + 2K) ,k,,%1%, + 


2k, kn Xp Hauck Dhak xa + 2h. k x)) + Top spec sa 


24°347273 24°44" 2*4 34°44%3 «1 
i. = Die oie wt a a 
2% 53 2 lis ee 22° 33%3 
o 2 “ e w e 
kak y + 2k 9% 1X5 + 2k 3% 1%3 -}- 2k 1% 1%4 + 2K 4X5%X3 + 


x. + ki ox. +k 


bed oe 
2k,,X,X, + 2k,,%_X, + 2x, (-k, + ky l 12*2 13*3 + 


24°24 34°37°4 
ky 4%,) + 2fago(% - x, (0)) + a4(Xq - ¥4(0))} (-k, + 


ki 0*1 + ky o%o + ko 3%3 + ko ,%,) + 2x) (-k, + k 3%] + 


KogXo tb Baatg © Kay%_) + Sie ee 


a3 (% - x4(0)) + an 4%,] (-k, b kia*1 + ko 4X + 
ey ee) 


The matrices in the first term are expanded and all terms are 


grouped in the form of (7.14). The resulting P coefficients are; 
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2 2 2 2 2 
eo RR eee © 995%7a t493%3q) + 


2 
WsX4q > + Ino*a%oq + 943*1a*3q + 914*1a%4a + 123*2G*3q + 
GonXoa%ua + 994*3a%4a t 2k, (65.5%5(0) + ay4X,(0)) + 2k, (a,9x,(0) 
+ 2 ,3%3(0)) 

2 


P, = 2k + 4b, Ky > Cdyy%yq + IyoXeg + 113%3q * 


tyy%ua? ~ yp (2o0%p(0) + 295%3(0)) = 2ky,(2,9%9(0) + 


a,3%3(0)) 


2 
Po m 72k, + Ab, kiko, = Cqyo%g + Ig9%oq + 993%3q + 


do4kqa) ~ pp (299%p(0) + a 3%3(0)) ~ 254,(a,,9%9(0) + 


a,4%3(0)) - 2k, - 2aok, ~ 2a,,k, 


Pa ~2k, + 4b k,k 


2 
3 4 *akaq ~ (4y3%1a t I23%2q F 9399%3a + 
934%4a) ~ 2 (299% (0) ¥ a95%3(0) = ky, (2,9%9(0) + 
243%360)] = eapagor SBA 3, 
¢ 

2 
By gt RAR So, “Meta © Uead.* Taare” 
Ws%4a? ” 294 (29%, (9) 4 23%3(0) . 2 yg (ayp%9(O) , 


a,3%3(0)) © 2k. ° 2a, ik, 


: 2 2 
Po = ¥ay * su - 2b, Ki, 
Pos = ka tilé.. -@b, 0° + We. + Ob be. + eek 
22 22 eee 4 24 12 22°22 42°24 
a? eee ee re 
33 33 5133 4 ~34 23°23 43°34 
a = £.. +14 - ope + oe, Ft Fae 
bly 4h sae 4 “44 34 44°44 
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Pp, =» Ik = que 


12 2 * yo ~ , *4"oq > Shi * 2559" ieee 

Pig 7 yg + ayy - 4, Kygkyy + 2aygky, + 2a,sky, 

Ti, xe Zk + ayy - 4bg Kygkyy + 2ky3 + 2e,aky, 

Pog 7 kag + dog ~ 4b, Kogkgy + 2Ky3 + 2aggkyy + 24k, 
Fe peter ta Bi3 tay 

Poy Ahoy + doy ~ 4b, Kogky, + 2ky, + 2kyy + agolo, 
+ 2a, ok, + 2a Koy 

3 ie 2s, + ayy > 4b, Kauk,, + 2g + 2aygky,+ 2a, sky, 
Ta 2a, ,k3, 


For each P coefficient to equal zero for ail values of the 
measured signals as required by (7.15), the following first order, 
linear differential equations must be satisfied; 

(B.5) 

. 2 2 2 2 2 

k = 2b,"k,” - Mayi*1q_ + 422%2a + 933%34 * G4 aa ? 

~ Wye*id*eq ~ 913%1d%3a 7 914*1a*4a ~ 9123*20*3¢ 
= do4%oa% na > 934% 3d%4a 7 22 (2q9%2(0) + ay 3%3(0)) 
- 2k, (a, 5%, (0) e a43%3(0) 

1 ot 2, "kjk, - 36411%14 + 412%24 + 113%34 * 914%4a 


: 19 (299%, (0) bs 253% (0) i Ky (249%50 us a,3%5(0)) 


~ e 
ii 


ye 
8 


2 
ee aa 3{41 2" 14 + do9%oq * 993%3q + I94%4a? 
“kyo (a,x, (0) + 24%, (0)) : k 9 (242% (0) ii a4 3%3(0) 


~ ky - agoks - ayoky, 
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/——e 
Kk, = 2b, kiko, - $43"14 + 4n3¥oq + 933%3q + 934%4q) 


- Kyg(2,9%(0) + a54x,(0)) - ks, (2,,5%5 (0) + a,4%,(0)) 
~ Aggky = aygky, 


k, : Zb, ky Kyy, z 5414" 1¢ + do4%oq * 194%3q + 944%4a 


= Ky,(255%,(0) + a95%3(0)) . Ky, 29% (0) + a,3%3(0)) 


- k, ” any ky, 
: 2 
kK * 2b, ky, - 4 q11 
— 2 
Koo = 2b, Koy, "4 Gog ~ 2Kyy > 28Q9Koo ~ 28, oko, 
2 
k33 = 2b, ka, "< 9133 ~ 289394 - 2845Ks, 
ee ee ~ 2k, - 2a,,k 
44 444 = q44 34 444d 
e e ] 
Kin = 2by Ky Ako, "3 dia 7 Kya 7 AeaK10 > B4oKi4 
ki3 ™ 2b, Kay Kay "4 di3 7 893ki9 > 843*1, 
k an 

1h Oy Rash ~* Set > Bis 7 Pye 
& 

Kos 2b, Koy ka, ‘ 2 493 7 ig > 8g9Ko3 - 83Ko07 Ay 0K 5, 

- ay gkoy, 
° 2 
sf a TS A "Ss dog ~ Kya 7 ¥o3 ~8aq%on 7 8nakug > 244*o4 
| a i Reta - : 934 7 33 7 823%, 7 843%yq > Aqa*s4 


where the boundary conditions at t = t, must be specified, 


£ 
Equation (7.3) is the expression for E at t = t-. Equating 


this to (7.12) evaluated at t = t,. gives 


£ 
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4 4 4 
(B.6) k(t) - 2 21 eee) xeGes) oh 2s Kap (tg) Malt pM (Eg) 


"5 { txlep + age pl" Fixe - x,(ep1} 


Expanding both sides and solving for the k(t ,)'s gives the boundary 


conditions on the k's: 
(B.7) k(t-) = 1 2 2 2 
£ y (f44%1q (tg) + FooX%oq (tg) + £55%34 (tg) + 
2 
fagkaq (Cg) 1] + fox g(t ed Xoglte) + fy yxy gltg)%gqlt_)t 
Eragkyqltpeyqlt_) + fogxoglte)xgqlt_) + 
foy%oq tg) eyqlt_) + f5,%3q (te) %4qCt,) 


k,(t.) = 
< (£1y% alte) + FyoXoglte) + Py5%3q) + Pry%yate?] 


nj 


k(t.) ) = 
2 iak [Ey o%qlt_) + Fyoxoglte) + Eogxgalte) + FopKy q(t] 


mie 


k(t.) = 
ae [Ey3%qlt_) + Fogxoglte) + Eggxgqlte) + Fay x,qlte)] 


mie 


k, Gt) a 
age [Ey palte) + fog%oqlte) + fay%qlte) + Fyy%yalte 


no} 


ki (t.) = 


nfm 
— 
fj 


Koo(te) = 


nfm 
ee] 
ih 


Kkag(te) = 


nope 
Ww 
vy 


Kyatts) = 


Nie 


kig(t,) = Lf, 


n | 


Kig(t,) = 


N}e 


kigft_e) = 


Nie 
Peis 
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